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Abstract 

In this paper we consider the dynamical system involved by the Ricci 
operator on the space of Kahler metrics. A. Nadel has defined an iter- 
ation scheme given by the Ricci operator for Fano manifold and asked 
whether it has some nontrivial periodic points. First, we prove that no 
such periodic points can exist. We define the inverse of the Ricci operator 
and consider the dynamical behaviour of its iterates for a Fano Kahler- 
Einstein manifold. In particular we show that the iterates do converge to 
the Kahler-Ricci soliton for toric manifolds. Finally, we define a finite di- 
mensional procedure to give an approximation of Kahler-Einstein metrics 
using this iterative procedure and apply it for blown up in 3 points. 



Introduction 

Let M be a compact Kahler manifold of complex dimension n. For any Kahler 
metric g we let uj := uig = ^f—Xjl-K ■ g^j{z)dz^ A dz^ denote its corresponding 
Kahler form, a closed positive (l,l)-form on M. For a Kahler form u; such that 
[uj] = ci{M) consider the space of strictly tj-plurisubharmonic potentials 

Aa[^] = {(y5 e C^{M) : uj + V^ddip > 0}, 

and the space ifeci of Kahler forms cohomologous to ci(M). 

In [Na], Nadel considers for a Fano manifold M, iterations on Ka^ defined 
inductively using the Ricci operator as follows. Let = Kuc^ denote the set 

of all smooth forms in ci(M) > and let Ric'-°^ be the identity operator. For 
any Kahler metric we let 

Ric(ci;) = ~ddlogdet{gij) 

denote the Ricci form of a;. It is well-defined globally and lies again in the 
ci(M) class. If it is positive we let Ric*-^-* (w) denote its Ricci form, and in a 
similar fashion we define higher powers of the operator as long as the positivity 
is preserved. The motivation for this construction comes from the simple fact 
that, when they exist, Kahler-Einstein metrics are by definition fixed points for 
this iteration process. Nadel asked whether these are all periodic points and 
proved the absence of periodic points of order 2 and 3. Furthermore he raised 
the question whether the existence of Kahler-Einstein metrics could be related 
to this iteration procedure. This question is also very natural. Actually, as we 
will explain later, one can define the operator Ric*^~^) using Calabi-Yau theorem 
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and see its iterations as a kind of naive discretization of the (normalised) Kahler- 
Ricci flow 

^ = -Ric{cut)+uJt (1) 
ot 

Thus, we expect some similar properties for the Kahler-Ricci flow and the iter- 
ations of the operator Ric'^"^^ 

We explain now the organization of this paper. We will show that some 
natural energy functionals arc decreasing along these iterations. This will give 
us a simple proof of the non existence of (non trivial) periodic points and thus 
answer Nadel's first question. Then we study the question of the existence of 
periodic points of infinite order. The behaviour of our dynamical system is 
closely tied with the existence of Kahler-Einstein metrics on the Fano manifold. 
We generalize Nadel's iteration scheme and define a family of natural operators 
RiCe, e > and see that their behaviour is particularly simple when the manifold 
is Einstein and e < 1. In Section 4, we investigate the case of Fano toric 
manifolds for which we can prove that the iterates of Ric*-"^^ do always converge 
towards a Kahler-Ricci soliton. This also gives a new proof of the existence of 
Kahler-Ricci solitons on Fano toric manifolds. We shall emphasise the fact that 
we don't use any flow or continuity method, contrarily to [Zlm2, WZ]. Then 
we relate the iteration procedure to the notion of canonically balanced metric 
studied by Donaldson in [Do.3]. This gives us an approximation procedure in a 
finite dimensional setup of the Kahler-Einstein metric when it does exist a priori 
on a Fano manifold. Wc apply our techniques to the case of blown up in 3 
points and give a numerical approximation of the Kahler-Einstein metric living 
on it. Finally we discuss the case of the manifolds with negative first Chern 
class and also non canonical classes. 



1 Positive and negative Ricci iterations 

Our first observation lies in the fact that Nadel's construction can actually be 
reversed. Let a be any form representing the ci(M) class. By the Calabi-Yau 
theorem [Yn] there exists a unique Kahler form in Kuc^ , which wc denote by lui, 
whose Ricci form equals a. Wc define the inverse Ricci operator by 

Ric<"^)(a) := LOi. 

Similarly we define higher powers Ric'"'^-' := Ric'"^-' o • • • o Ric^~"'^\ making 
repeated use of the Calabi-Yau theorem. This leads us to consider for a given 
Kahler form ujq the sequence 

Ric(ti>j+i) — LOj (2) 
for an integer j > 1. This can be also written 

^■^ - = -RicK+i) 

which justifies naively why our iterative procedure (2) should have some simi- 
larities with the Kahler-Ricci flow (1). 
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For a positive k we will denote Ka^i the maximal domain of definition of 
Ric'-''^ Wc obtain naturally a filtration of Ka^ci ■ We also let Ka^ci denote the 
set of all metrics in Kuci whose image under Ric^ ' is nonnegative. 

The advantage of this construction lies in the fact that iterations improve 
the positivity of a. Indeed, if we consider the filtration {Kai{^}j>Q to measure 
positivity in some sense, then the image of Ric'"''^ hes in Ka'^l^ and we may 
iterate the inverse Ricci operator to any desired power. Note that this yields 
homeomorphisms 

Ric(^) : (Aa«, II • W^^.,) ^ (Aa^'+^'Ml • 11^™-..) 

for all (3 e (0, 1), m, ^ e N U {0}, k € Z such that I + k > and m > 2k. We 
let Kai^' denote the set of all L°°(M) limits limfc_oo Ric^'^'^a, for a € ife^^ 
(when they exist). 

Finally, let G be any connected compact subgroup of the group Aut(M, J) 
of holomorphic diffcomorphisms of (A/, J). Denote by Kan{G) the space of G- 
invariant Kahler forms in Kan- Such forms exists as can be seen by averaging 
over orbits of G with respect to the Haar measure of G. We remark that Ric'-"^-' 
maps Kaci (G) i^ito itself. 



2 Energy functionals on the space of Kahler po- 
tentials 

In the Sections 2.1 and 2.2, we present some properties of some well-known 
energy functionals. Let Vl E H^{M, M) denote a Kahler class. We call a function 
A : Kan x Kan — > Hi an energy functional if it is zero on the diagonal. By an 
exact energy functional we will mean one which satisfies in addition the cocycle 
condition (see [Ma]) 

A{uJi,U2) + A{LL>2,t^3) = A{uJi,UJ3). 

In all the paper, V will denote the volume of the manifold with respect to 
[w], i.e V ~ Jj^j The energy functionals /, J, introduced by Aubin in [Aiif], 
are defined for each pair {uj, uj^p uj -\- \/—lddip) by 

n— 1 

V Jm y Jm 

J{cu, cu^) = \ f V^d^ A A Y.{n - i)u:^ A (4) 

We note some of their basic properties for which we refer the reader to [Aul, 
Si, Til]. Note that /, J and I — J arc all nonnegative and equivalent. One 
may also define them via a variational formula. Connect each pair [oj^lo^-^ :— 
uj + \/ —\ddtpi) with a piecewise smooth path {wi^t}. Then, e.g. for I ^ J, we 
have for any such path 

{I- J)[u,LOp,) = -^ j ^Atipt^^Adt. (5) 

y J[0,l]xM 
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2.1 The functional 

Let 

F^{LOo,^) = -{I-J)i^o,Lo^)-^f V^l- (6) 

For a Kahlcr manifold of positive or negative Chern class define the following 
functional on Aa^d x Ka^a 



M 



with fj, = ±1 respectively. For a Kahler manifold of zero Chern class define on 
Kan X Kan the functional 

^ JM 

The critical points of the functionals are the Kahler-Einstein metrics. They 
are absolute minima [Di, Til]. Indeed the second variation of F at a critical 
point in the direction of the plane spanned by t/"!, ■02 S T^Ka^^^-^ is given by 

^ / [^5^(VV'i,V02)-mV'iV'2]<, M = ±1,0. (7) 

V J M ^ 

This is seen to be a strictly positive (0, 2)-tensor on Ka\^^^ for ^ < —1. For 
/Lt = 1 it is nonnegative and vanishes precisely when tpi and 'ip2 are proportional 
and eigenfunctions of eigenvalue —1 of Ag (see [Til, p. 64]). In that case this in- 
finitesimal variation corresponds to holomorphic automorphisms and to moving 
within the set of Kahler-Einstein forms. 



2.2 K-energy and Ek functionals 

The Chen-Tian energy functionals Ek, k — 0, . . . ,n, are defined in a similar 
manner by 

Ffc(c^,c^^J = / A^,v^tRicKJ'=Ac.;-'=Adt 

^ J[0,l]xM 

_ / ^*(RicKJ'=+i - (8) 
" H"([M]) ■ 

This gives rise to well-defined exact energy functionals independent of the choice 
of path [C'T]. The K-energy, Eq, was introduced by Mabuchi [Ma]. The follow- 
ing formula is taken from [TiO, Section 7.2]. For the derivation (for any Kahler 
class) we refer to [Ch] where an equivalent expression is given. 



Proposition 2.1. Let h be a function satisfying Ricw — jiuj = y/—lddh. One 
has 

1 /■ 1 /■ 

Eo{u;,io^)^- log-fy;^^iI-J)iu;,Lo^) + - ~ lu^). 

^ JM ^ ^ JM 
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Remark that the Ek functionals vanish on pairs joined by a one parameter 
subgroup of automorphisms through the identity [CT, Corollary 5.5]. 

Definition 2.1. We say that an exact functional T is bounded from below if 
T{u>, ujip) > C for every ui^ £ Kan. We say it is proper on Kau{G) (in the sense 
of Tian) if there exists a function p : M ^ M satisfying \mis—,oo p{s) — oo such 
that T{bj,bj^) > p{{I - J){LO,uJtf,)) for every lo^ € Kan{G). 

This is well-defined, in other words depends only on [w] since the failure of 
/ — J to satisfy the cocycle condition is under control with respect to the two 
base metrics, 

[I - J)[UJ,UJ^,) ~ {I - J){W^,,UJ^,) = {I - J)[UJ,UJ^,) ~ y / (^1 (w;^, -UJ"^,). 

^ J M 

2.3 A lower bound for the energy functionals 

We now study the iterations of the Ric^^ operator for Fano manifolds by 
analysing the behaviour of the functionals that we have introduced. Firstly, 
we recall a result of non-negativity of Bando and Mabuchi, 

Theorem 2.1. [BM, Theorem A],[Ba, Theorem 1],[SW1, Theorem 1.2]. Let 
{M^loke) be a Fano Kdhler- Einstein manifold. Then for i = 0, 1, 

for all CO € Ka^ with equality if and only if lu is Kdhler- Einstein. In that case 
there exists a holomorphic automorphism homotopic to the identity h such that 

h^LUKE = 

We give a sketch of the proof for i — with an emphasis on the features 
that will be useful in later sections. Consider the deformation {tiJipt} C Ka^ 
constructed from two paths, solutions of the following Monge- Ampere equations 

< = e*^+^'c^",t e [0,1] (10) 
= e^-(*-i)«w",t e [1,2] (11) 

where Ricw — u = \/—lddf with the normalisations 

/ e*^+"'tj" = / ef-'^'-^^'f'Lo'' ^ V. 
Jm Jm 

Note that the first path is the one used in Yau's continuity method proof [ V;i]. It 
connects any point lu in Koc^ to Ric''~^'cj in Ka^^^K The second path, introduced 
by Aubin in [Aul], is used to connect any point in Koci to a Kahler-Einstein 
metric. 

The existence of the first path is equivalent to the Calabi-Yau theorem. The 
second path may not exist in the presence of nontrivial holomorphic vector 
fields but Bando and Mabuchi show that arbitrarily close to uj in the C°°{M)- 
topology there exist metrics for which such a path exists. Since the K-energy is 
continuous this will be sufficient for the argument (Cf. [BM],[SW1, Section 3]). 

Now, for t G [0,1] one has 

Ricw^ = (1 - t)Rictj + tw, (12) 
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and 
hence 



At^ = / + ct, (13) 



d 1 /" 

—Eo{ujKE,i^.pt) = -y J y^CR-^cuj 



•fit - ^vt) A nuj^ 1 



1 



V 



M 



(14) 

with t £ [0, 1], from which wc conclude 

Eii[uJKE,l^^t) < E(3{uJKE,t^)- (15) 

Next, for [1,2] 



and 
hence 



RictJ^ = (2 - t)Lj + (t - l)w^, 



d „ , . 1 



^-EoiujKE,i^<Pt) = J "^(-(2 - t)v-199w^) A ntj, 



n-l 



-(2 -<)-(/- J)(^K£;,w«) <0 



(2-<)-i / ((A,^)2+t|9<^|2)a;;<0, (16) 



'M 

where t G [1,2]. The theorem now foUows in this case. Song and Weinkove 
extended this argument to Ei using two detailed computations. The first shows 
that while Ek may not necessarily be monotone (when the path exists), one still 
has Ek{ujKE,^ipi) > Ek{ujKE,^,p2) = 0- other words, 

Theorem 2.2. [SWl, Theorem 1.1] Let {M,ujke) be a Fano Kdhler-Einstein 
manifold. Then for any lo G Kal?^ and for each k = 0, . . . ,n one has 

Ek{uJKE,tL!) > 0, 

with equality if and only if to is Kdhler-Einstein and h*LOKE = ^ with h a 
biholomorphism homotopic to the identity. 

The second calculation shows that when fc = 1 , one has 

Ei{uj^^,uj) > Ei{uJo,uj) = 0. 
Explicitly, their computation shows that 

Ekiuj^, ,uj)= y Jj^j V^difii A Blpi a '£'^=0 ^tUj' A 



h Im Eti i:X\)fi^9dfy A c"-, (17) 



V JM Z^i=l \k+ 

with a, = '""''1^,'+^^ if < i < fc - 1 and ai {i k < i < n. Since the 

' n+l — — ^ n+\ — — 

last term is positive on Kac^ for fc = 1 they conclude their proof. 
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2.4 A system of Monge-Ampere equations 

Now, let Lu = luq denote an initial Kahler metric for our iterations. We present 
the iterative procedure defined by (2) in terms of Monge-Ampere equations. Let 
ifi be a Kahler potential with 

Ric(tJo + V—iddtpi) = ujo. 

Set the Ricci deviation h h^^^ of luq as \/^lddh ~ Ricwo — wo- The function 
h thus given is for the moment determined only up to a an additive constant. 
The equation then becomes 

— -s/^SS log det g^^ = ojo = — \/— Ti991ogdet g — ^/~^ddh (18) 



or 



that is 



-199 log 



-Iddh 



n h n 

w^pi = e w 
together with the volume normalisation 

This determines ipi only up to a constant, which will be fixed in the the next 
step. Put uji ~ Lu^-^ . In the second step we solve 

Ric([Ji + V^99(^2) = 
and uji — Ricwi = lui — ujq = \/ —Iddipi. The Monge-Ampere equation is now 

with (fii determined uniquely by 

/ e'^^'^^w" = 1. 

Iterating this procedure we have Ric^~'''a; = uj^i for each I £ N where 



and each of the tpj is uniquely determined by 

— / e''"^j=i'^^w" = 1. 



From now on we set 
and 



uji = cj^j = ojo + V—Tdd^i 



(19) 



(20) 
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2.5 Monotonicity of the energy functionals 

The following proposition describes the monotonicity of the K-energy and 
along the iteration. Note that the equivalent for the Kahler-Ricci flow is well- 
known. 

Proposition 2.2. Let /i = 1. Then 

^ Jm ^ J M 

Ei{io,uji) < 0, 

with equality if and only if lu is Kdhler- Einstein. Furthermore all the Ek with 
k = 2, . . . ,n decrease along the iteration starting from the second iteration. 

Proof. To show the first inequality we note that 

En{ujk-i,(^k) = ^ 1 -ipk~iuJk - {I - J){uJk-i,uJk) - ^ 1 ipk-i{uJk-i - ^k) 
= -{I - J) {iJk-1, ujk)- ^ [ Vk-i^l-x- (21) 



V 



M 



The first term is nonpositive with equality if and only if uj^ — ujk~i = Ricwfc, 
while the second term is nonpositive since 

1 ^ / = ^ / ^^'"'-'^k-i >Yr I (1 - ^fc-i)c.Li. (22) 
^ Jm ^ JM ^ JM 

Since 

I 

Eo{uj,uJi) ^ Eo{ujk-i,uJk), 
fc=i 

the conclusion follows. 

The second inequality follows similarly, indeed 

fi(tjfc_i,Wfc) = -(/- J)(wfc-i,Wfc) - / (fiki^k^ 

^ JM 

and with (22). The third inequality follows from 

Ei{ujk~uuJk) = 2Fi{uJk-i,ujk) + ^ / ipk{i^k + ojk~^ ^ ujk-i) 

^ JM 

from the formula relating Fi and Ek. Since both summands are negative, we 
conclude (note that uik-i — u>k = —V—^ddipk). Finally, the decrease of the 
functionals are proved using (15) and (17). □ 

Remark 2.1. For the rest of the energy functionals, we get by direct computa- 
tions that with respect to a Kdhler- Einstein metric uj, 

Ekito^LO^) = {k + l)Fi{Lo,LO^)-^ f f^^{ul + ...+ujl-'' hiRicu^f) 

V JM 

= {k + \)Fi{uj,uj^p) - Bk{uj,Lo^p). 
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where Ricujip = uj^ + \/—lddfui ■ Note that when the Ricci curvature is positive 



we have 



V .Im 



1 

'v 



M 



+ ^ I f^^u;^Aiuj:;-' + ...+u;;-'^A{Rtcu;^r~') 
< Bo(w,cj^)<0. (23) 

3 The dynamics of the Ricci operator 

We are now ready to answer the question raised by Nadel [Na] . Note that some 
of the results presented in that section were discussed in details by the author 
and Y. Rubinstein [Ru] during the period^ of their collaboration. 

Theorem 1. Let (M,uj) be a Kdhler manifold with positive first Chern class 
and assume that Ric('')(w) =uj for some fc e N. Then cu is Kdhler- Einstein. 

Proof. Note that the nonexistence of fixed points of negative order implies that 
of positive order, and reversely. Therefore assume that for some uj £ Ka^ and 
some / S N one has Ric^~'^((jj) = uj. By the cocyclc condition wc therefore have 

(-1 

= £;o(u;,Ric(-')a;) ^ £;o(Ric(^*)a;, Ric^^'^^^w). (24) 

1=0 

On the other hand, from the first part of (11) 

So(c^,Ric(-i'c^) = -i / (l-t)(A^^^)2c^" Adt-(/- J)(c^,Ric(-i)cc>) 



V 



A/x[0,l] 



Thus Eo{uj, Ric'"^-*!!;) < 0, with equality if and only if Ric'^^'tj = cj. Therefore 
each of the terms in (24) must vanish identically and we conclude that (Af , uj) 
is Kahler-Einstcin. □ 

Proposition 3.1. Let M he as above and assume that uj G Kac'f^ for all k N 
and that uj is not Kdhler-Einstein. Then limfc_,oo Ric'^^'w does not exist in Ka^ci ■ 

Proof. If ujoc = limi^oo Ric^'^w exists and is smooth it satisfies Ricwoo = ^oo- 
But £'o(ijJ, ijJoo) > contradicting (2.1). □ 

Let G denote the Green function for A = Ag with respect to (A/, w) with 
G'(x, y)uj"{y) ~ and A{uj) = — inf G such that 

/(x) - i / /c." = -i / G(.T,2;)A/(y)c."(y), V/eC°°(Af). 

Then, one has the following estimate due to Bando and Mabuchi. 



1 September 2003 to December 2006. 
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Theorem 3.1. [DM] One has 

A{uj) < ydiam(Af,a;)2. 

//Ric(cL;) > ecu for some e > then diam(M, a;)^ < ^ by Myers' theorem. 

As an immediate corollary we have 

Lemma 3.1. Let M be a Fano manifold. Assume that the K-energy is proper. 
Let (ci;;);gN be a sequence of Kdhler forms on which the K-energy is bounded from 
above and such that there exists Io{ijJq) G N and e > with Rica;; > eLUi,yi > Iq. 
Then there exists a constant Ci depending only on (M, wq) such that 

ll<i'/||L-(Af,^) <C^i, VZeN 

where lui ~ luq + \/—ldd^i and y Jj^^ e''"*' ^ = 1. 

Proof. Let Gi be the Green function for A; ~ ^8 uji (i-*^-; the Laplacian with 
respect to {M,uji)) satisfying ^^jGi{x,y)ujl{y) = 0. Set Ai = -infMxMGj. 
Since — n < Ai$/ and n > A/$/ the Green formula gives 

Mx)-^f^iu;^^ -i/,,Gi(x,2/)A<i>;(yK(y)<nAi, (25) 

-yJ -^i^'i - j^,Gi{x,y)^<i>i{y)u;l{y) > -uAi. (26) 

Hence 

oscM^i<n{Ai+ Ai) + I{uj,uji). (27) 

Since Eq is proper on Kac^{G) in the sense of Tian, if Eq{loq^ ■) is uniformly 
bounded from above on a subset of Kaci{G) so is I{uj,-). We conclude that 
I{uJo,U!i) is uniformly bounded independently of I. Finally, the Proposition 
follows from Theorem 3.1 which provides a uniform bound for Ai. □ 

As a consequence of the properness of the K-energy for Fano Einstein man- 
ifolds [Til, PSSW] we obtain the following, 

Corollary 3.1. Let M be a Fano Einstein manifold with no nontrivial holo- 
morphic vector field. Consider the sequence of Kdhler metrics (ci;;);gN defined by 
the system of Monge- Ampere equations (19) and assume that for I sufficiently 
large there exists a constant e > with 

Ric{uJi) > eujQ. 

Then loi converges to the Kdhler- Einstein metric when I tends to infinity. 

Finally, inspired by [Pa], we derive the following. 

Proposition 3.2. Let M be a Fano manifold with G a maximal compact sub- 
group o/ Aut(A/, J). Consider the sequence of G -invariant Kdhler metrics LOi 
defined by the system of Monge- Ampere equations (19). Assume that there exists 
a constant 1 > k > such that for I sufficiently large, 

(2 - k)oj''(, > LO'^ > KUjj; 

where loq is G-invariant Kdhler metric. Then AI is Kdhler- Einstein and oji 
converges to a G-invariant Kdhler- Einstein metric when I tends to infinity. 
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Proof. Thanks to the proof of Lemma 3.1 and Theorem 3.1, we are reduced to 
prove an upper bound for I{uj, tui). But if we denote ^'1' (x) = sup{0, and 
= mf{0, ^i{x)}, we obtain 



^ "'{w">Wo} ^ J{w,',''>w,"} 

< (1-k)0SCa/$; 

Together with 

oscM$i < n{Ai + Ai) + I{lu,lui). 

this gives us to the C° bound for . Now it is a standard argument of Monge- 
Ampere equations to get the convergence in the C°° topology. □ 

For a Kahler manifold M, we consider the family of Monge- Ampere equations- 
(11). We introduce the Aubin operators Ric^ by setting 

RiCe(tj) UJ,fi + , 

for each e G [0, 1] such that (fi+e, solution of (11), exists. Note that Rico(a') = 
Ric^~^-'tj and Rici(a') — ujke- Formally, one can think 

RiCe = ^^-^(Ric - eld) 

and that we have defined the following sequence of Monge- Ampere equations 
(tJo + V^dB^jY = e(^-i)*^-i-<^*^ {uja + \/^o>9$j_i)" (28) 

Let G C Aut(M, J) be as before. Then from uniqueness of solutions of the family 
of the Monge- Ampere equations (11) we conclude that RiCg maps Kaci{G) into 
itself. 

We recall the definition of a-invariant introduced by Tian. 



sup{a > : sup / e-°"^uj" < oo}. (29) 



In [TiO] it is proved that aM is a positive holomorphic invariant of Fano mani- 
folds. Regarding the existence of the Aubin operators we state the following 

Proposition 3.3. 

(i) [TiO] Assume that M is Fano. Then the operators RiCe exist for all e G 
[0,min{l,S±iaM}). 

(ii) [BM, Theorem 5.7] Assume in addition that the K-energy is bounded from 
below. Then the operators RiCg exist for any e G [0, 1). 

We now recover by a conceptually simpler method a theorem of Tian [Tit]. 

Corollary 3.2. Let M be a Fano manifold. Let G be a maximal compact 
subgroup of Aut{M, J) and assume that the K-energy is proper on Koc^iG) and 
let e G (0,1). Then there exist G-invariant Kdhler- Einstein metrics. All such 
metrics are the limit points of the iterates o/ Ric^ on Kac-^{G) in the C^{M)- 
topology. 
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Proof. Since the K-energy is proper on Kuc^ (G) and in particular bounded from 
below, RiCe are defined for each e G (0, 1). As before we obtain (27) with cui = 
Ric^'^(a;) with each cui in Kac^iG). Since by (15) and (16) the K-energy decreases 
along iterates we still have a uniform bound on I along the orbits. By Theorem 
3.1 we also have a uniform bound (depending on e) on AiJMc^P [uj)) . For each 
/3 e (0,1) the C2'^(M)-estimates now follow from [Au3, Chapter 7] and the 
higher order derivatives follow from these, by bootstrapping. We may therefore 
extract a converging subsequence Ric^'^''(u;) in the C^'''(M, w)-topology whose 
limit Woo lies in C^(M). Moreover, since the K-energy is bounded from below 
we have 




As each of the summands is nonpositive one has i?o(woo, RiCe(a;oo)) = and it 
follows that Woo G KuciiG) is smooth and Kahler-Einstein by (15) and (16). 
Since this is true for each converging subsequence we conclude that, in fact, the 
sequence of iterates itself converges. □ 



In fact we also get, 

Theorem 2. Let M be a Kahler-Einstein manifold with positive first Chern 
class. Let G be a maximal compact subgroup o/Aut(Af, J). 

Then for any e G (0,1), uj e Kac^{G) there exists a biholomorphism h of 
{M,J) such that one has limj^oo R-iCe'^'i' ~ h*ujKE in the CP" {M) -topology for 
some G-invariant Kahler-Einstein form loke- 

4 The case of toric Fano manifolds 

On a toric Kahler manifold M with positive first Chern class, a result of X-J. 
Wang and X. Zhu [WZ] asserts that there does always exist a Kahler-Ricci soli- 
ton, which is unique up to holomorphic automorphims. A Kahler-Ricci soliton 
is a pair {X,LOg) where X is a holomorphic vector field on M and uj = ojg & 
Kahler form that satisfies 

Ric(w) = w-hLx(w) (30) 

where Lx is the Lie derivative along X. In fact, the Kahler-Ricci soliton in 
the toric case is a genuine Kahler-Einstein metric if and only if the Futaki 
invariant vanishes which implies that the holomorphic automorphism group of 
the manifold is reductive. 

Given a nontrivial holomorphic vector field X on M, it is well known by 
Hodge theory that there exists a unique smooth complex-valued function 6x 
such that 

ix{uj)^ v^BOx 

with J^e^^LO^ = V. Let h S C°°[M) be the Ricci deviation of uo determined 
by Y Im s'*'^" — 1- From [WZ, Lemma 2.1], there exists a unique holomorphic 
vector field X such that 

Fx{v) := / v{h - 6lx)e''^w" 
Jm 
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for all holomorphic vector field v. Note that Fx is actually independent of the 
choice of uj. For such a vector field X, to find a Kahler-Ricci soliton {X,lu + 
\/~\ddip) satisfying (30) is equivalent to solve the Monge- Ampere equation 

In [Zhul], an existence result (proved by a continuity method argument) is given 
for the Monge- Ampere equations of the form 

det(5,- + v^,-)=e/-*«--*^(-)). 

where / is smooth and t varies from to 1 (note that t = corresponds exactly 
to Calabi-Yau theorem). This can be seen as consequence of the work of S. 
Kolodziej [Ko] for solving Monge- Ampere equations, and the following result 
[Zhul, Corollary 5.3]. 

Lemma 4.1. Let (AI,lu) a Kdhler manifold with a non-trivial holomorphic vec- 
tor field X . Suppose that ip is smooth, u -\- y/—ldd(p > and X(ip) is a real- 
valued function. Then there is a uniform constant C independent of such that 

I^MI < c. 

Hence, we can define naturally the sequence of Monge- Ampere equations 

which turns out to be equivalent to solve for all j > 1 

Ric(tjj) = tjj_i + Lxi^i) (31) 
for LOj = ti)j_i + \J~^ddpj =0^0 + V— 199$^ . 



The goal of this section is to prove that under this setting, one can always get 
a C*^ estimate for the potentials <f>j, which means that the iteration procedure 
defined by (31) converges to the unique Kahler-Ricci soliton on the manifold. 
We follow the techniques developed in [\A Z] using real Monge- Ampere equations 
to derive the uniform estimate for <I>j. 

We shall need the following lemma that we will apply with the convex poly- 
hedron f2* associated to our toric manifold. 

Lemma 4.2. LeA, a bounded convex domain in M". Then there is a unique 
ellipsoid £ called the minimum ellipsoid of O which attains minimum volume 
among all ellipsoids containing fi, such that 

-Ednas 

n 

where r£ denotes the dilatation of factor r of £ with concentrated factor. 
Note that the vector field X can be expressed in the form 

X = Y, c,Xi 

i 
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with Xi = for = Xi + \/^l^i)i=i,--n afBne logarithm coordinates on 
the maximal torus T of Aut(M, J). In [WZ, Section 2] it is shown that these 
constants Ci satisfy the relations 

/ y^e^'=^'"y'dy = (32) 
Jn* 

for all i — 1, ..,n.. Furthermore, there exists a Kahler form ujq, invariant under 
the maximal compact abelian subgroup A'o of T, and determined by a convex 
function wq on R", i.e on T 

LUo = ^/^ddwo- (33) 
In all the following, we shall denote det(wo) := det ((wo)y")- We obtain 

lUq ^ — dct{wo)dxi A ... A dxn A d9i A ... A d0n. 

Note that the function wq is determined uniquely by the convex polyhedron 
n*. Indeed, if {p^^''}i=i,..,m denote the vertices of fl*, wq is the convex function 
defined by wq ~ log (X^IILi e^^'"^^) . After normalisation and without loss of gen- 
erality, we can assume dct(wo) = e~^~'^° . We initialize our iteration procedure 
by chosing this way the Kahler form wq. Let's call Wj = ifo + $j with a 
^Co-invariant potential. Now, there exists a constant cx such that 

ex+X(<i>,) = ^c,^ + cx. 

This means that to solve (31) is equivalent to solve the real Monge-Ampere 
equation on M", 

det(wo + $j) = dct(wj-) = e""-^"""^-'"^"^ (34) 
Lemma 4.3. Define 

rrij ~ mi Wj . 

Then there exists a constant C independent of j such that 

mj < C. 

Proof. For an integer fe, we choose Xj^k a point in R" such that Wj{xj^k) = nij+k. 
We introduce for j > 0, the set 

Akj = {x e R" : mj + k < Wj{x) < m-j + k + 1, 

rrij+i + rj+i,k < Wj+i{x) < nij+i + rj+i^u + 1} 

where r^+Lfe = [wj+i{xj^k) - "ij+i] + k. For any integer > 0, Uo<j<j,yljj 
is convex. On the other hand, we know that _Dwj(R") = Vl* and thus Akj is 
bounded. Let's define dmax = supj^g^j. {c;j/;}. We get on Ao.j, 

From Lemma 4.2, there exists a linear transformation y = Tr{x) which leaves the 
center of the minimum ellipsoid of Ao.j invariant such that i?j^/„ C Tr(Aoj) C 
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Bji for a certain well-chosen constant R. Let's try to obtain some information 
about this radius. Define 

Vj = \ (e-^--''— " [\y - - ) j + + '^j+i.o + 1 

where Dj the center of the mininmm ellipsoid of Aqj-. Then, in Tr(ylo.j), 

det(wj) = e-^:^-'^— ig-m, 

and of course, > Wj+i on (9rr(Ao.j). The comparison principle for Monge- 
Ampere operator gives that Vj > Wj+i on T'r(Ao.j). Hence, 

m,j+i+rj+i.o < w-i+i{y-i) 



This gives us the following inequality with Co = \/2n (^e~^-'^~''"'"'~^) ^" 

i? < Coe^ 

On the other hand, since Tr{Ak,j) C -B2(fc+i)_Ri we obtain 

k 

< Vol{S"-^)'^e-"'^-''2{k + l)R" 

k 

But, on another hand, we also have thanks to (34) and by transformation y — 

Dwj+i{x), 

[ e-"'^ da; = e'^^ ( e^' dy = Cq. (35) 
Jr" Jn* 

This gives the expected inequality. □ 

Lemma 4.4. Let Xj^ be the minimum point of Wj . Then there exists C > 
such that \xj^Q\ < C for all j . 

Proof. The fact that < sup^jg^j. \x\, and the previous lemma gives the ex- 

istence of R' independent of j such that inf 93^,(2.^. ^-j Wj < nij + 1. By convexity, 

\Dwj{x)\>l/R' (36) 

in R^\BRr{xjfi). Hence with (35), for any e > 0, there exists Re sufficiently 
large such that 



/ 



e-^'^dx < Ci 



Jx-Xj^o\/R' 
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where i?c is independent of j . 

Now, assume that \xjfi\ is not bounded. For any e > 0, there exists a large 
constant C > such that if \xj,o\ > C then, 



> - inf \x\ 

2 xedn* 



for X e Bfi^(xjfi) and where we have set ^ 



This can be viewed by 



considering the restriction of Wj+i on the ray Oxj^ and by its convexity together 
with the fact that Dwj+i{W) ^ n* 3 {0}. Now with (35), 



Ib, 



Furthermore, we have 



^e-^>^ inf 



(37) 



dwj+i 



9C 



< 



sup \x\ 

xen* 



"■\BR,{xj,a) 



With (37) and (38), one obtains for e sufficiently small that 



Now on another hand by (34) and (32), 



>dx > 0. 



< e sup \x\ (38) 



(39) 



dwj+ij2 



dx. 



dxi 



"'dx 



which leads us to Jjg„ ^^q^^ e — and hence a contradiction with (39). Thus, 
\xjfi\ has to be bounded. 



□ 



Lemma 4.5. One has the upper bound 



sup(f>j < C" 



M 

for a constant C independent of j . 

Proof. As explained in the proof of [\A Z, Lemma 3.4], by convexity, it is suf- 
ficient to give an upper bound of <i>j(0). Since \xjfl\ < C by Lemma 4.4 and 
\Dwj\ < sup^gQ. |a::|, we obtain Wj{0) < C . With (35), we get |wj(0)| bounded 
and hence |<J>j(0)| bounded. □ 

Lemma 4.6. One has the lower bound 

inf $, > C" 



for a constant C" independent of j . 
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Proof. Let {p^'^}i=i,..,m denote the vertices of il* and define 

v{x) — max {{x,p^^^)}. 

i—l...,m 

Then the graph of u is a convex cone with vertex at the origin and is an asymp- 
totical cone of the graph of wq, i.c one can check that the foUowing incquahty 
is satisfied 

\v - Wo\ < Co. 

Let's denote 

Zj{r) = sup{v-Wj){x) 

\x\—r 

and assume that the supremum is attamcd at Pr.j- Then with ^ = y^^, we get 

z'^{r)<d^{v~w,){pr^,). (40) 

Claim. There exists tq independent of j such that for r > ro, Zj{r) < Ar^^ for 
a constant A independent of j. 

If the claim is proved then by integration, supj^ggn {v — Wj) is bounded inde- 
pendently of j and the lemma is proved. Let's now prove the claim. Let F^'^^ 
be the faces of the graph of v and suppose that 

F^^^ d {xi > c'o\x\} n {Xn+l ^ 0} (41) 

which is possible by a change of variables (cg is here a well-chosen positive 
constant and x = (x2, a;„)). As explained in [WZ, Lemma 3.5], if the claim 
does not hold, there would exist a point qrj in {xi = ■^Prj} such that 

Wj{qr.j) = mf{wj{x) : x G F''^^ n {xi = —Pr.j}}- 
It follows from (40) that for a certain constant Cg, 

Wj{qr^j)-Wjipr,j)>^. (42) 

Define for a convex function / the set 

Nf{x) = {p e R" : f{y) > f{x) + p ■ {y - x) Vy e R"} 

and for a set S, Nf{S) — UxesNf{x). Then, for ip defined on R" such that 
its graph is a convex cone with ^PiPrj) = ^jiPr.j) and ip ~ Wj on {x G R" : 
Wj{x) — Wj{qr.j)}, one has 

N^{S„) C N^^iSo) (43) 
where we have defined So = {x E R" : Wj{x) < ■Wj{qrj)}. From (42), one gets 

dist(0, So) > C2r 

and by (36) applied at step j — 1, we obtain for x Sq the existence of a constant 
C3 > such that 

Wj^i{x) > C3\x\ — C > C2r — C. 
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But now, with (34), there exists a positive constant Cg such that 



Vol{N.^^{So)) < f detK) < C3yo/(5o)e-^3- 



J So 



and on another hand, with (41) and (42), there exists C4 such that 



Vol{N^,{So)) > 



,n+l 



which leads us to the expected contradiction because of the inclusion (43). □ 

Once the C° estimate in hand for $j , one can derive C°° estimates by using 
the results of [Ya , TZ]. Finally the results of this section gives the 

Theorem 3. Let M be a Fano toric manifold with Kdhler-Ricci soliton {X, Woo)- 
The iterations defined by 



for i > 1 and with ujq given by (33), converge to a Kdhler-Ricci soliton {X,uJao). 
If the Futaki invariant vanishes, then the iterations converge to the Kdhler- 
Einstein metric on M . 

Remark 4.1. In fact our theorem gives another proof of the existence of Kdhler- 
Ricci solitons on Fano toric manifolds without using any flow or continuity 
method. 

5 Applications and numerical results 
5.1 The finite dimensional picture 

In [Do3] , Donaldson has introduced the notion of i/-balanced metric for a fixed 
volume form v and proved its existence [Do3, Prop. 4] under some very general 
conditions [Do3, p. 10]. These metrics have the properties to solve the Calabi 
problem, i.e to converge towards the Kahler metric that has volume form v in 
a given Kahler class. 

Let us fix a volume form v on a smooth projective manifold M with a 
polarisation L, and choose r G N sufficiently large such that M is embedded by 
the holomorphic sections of in the projective space PiJ°(M, L^). We set iV^ — 
h^{M, U) which is finite since M is compact. One defines a j/-balanced metric at 
rank r as the fixed point of the map : Met{H°{M, U)) Met{H°{M, U)), 



where G is a hermitian metric of H'^{M, U') and (S'i)i=i,..,7Vr is an orthonormal 
basis of II^{M, L^) with respect to G. Donaldson proved that the compositions 
of the map give a convergent sequence of metrics in Met{H^{M,U)) (and 
thus on Met{L^) by the Fubini map FS, see [Do3, p.4, Section 2.2]). The limit 
is called the j/-balanced metric. 

Notation 5.1. For a smooth hermitian metric h € Met{L) on the line bundle 
L, we denote Ci{h) G 27r[L] its curvature. 



Ric{ujj) - Lx{i^j) 



'3-1 
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Theorem 4. Under these settings, let's call G AIet{H^{M, U")) the sequence 
of v-balanced metrics of order r. Then ci{FS{HrY/^) converge to a Kdhler form 
Woo in [ci(L)] that satisfies 

Proof. To prove this theorem, we use the powerful Calabi-Yau theorem. Hence 
we know the existence of a Kahler form lo in [ci(L)] such that = v. We use 
X. Wang's theorem [Wa] with the trivial bundle and L. There is a Hermitian- 
Einstein metric on these bundles and the metrics Hr are 'balanced' with respect 
to LO in the sense studied by Wang. This is due to the obvious fact that the 
considered bundles are Gieseker stable. Thus, one obtains directly the conver- 
gence of the sequence of metrics FS{HrY/^' <E Met{L) to the metric h^ with 

Cl{hL)=LU. □ 

We now assume in this section that M is Fano and consider the polarisation 
L = —Km- We fix an integer r sufficiently large such that M is embedded by 
the holomorphic sections of L'' . Let's consider a smooth hermitian metric ho on 
L with ci(/io) = ujQ and let's call f^jg the Ricci deviation of ujq. Now for each fc, 
one can define a Ric'-"'^-' (u;o)"-balanced metric at rank r in the following way. 

Consider Hilb^glh^) the L^-metric induced on the space of holomorphic 
sections H^{M^L'^) with respect to wq, given by 

Hilb^,{{,)){S,,S,)^ [ {S,,S,)ef-ouj^. 
Jm 

Now, for a given hermitian metric Hq on H^{M,L^), we define the metric 
FS{Ho) on U by 

where the (S'i)i=i...,jVr form an _ffo-orthonormal basis of H^{M, L"^). Then, from 
[Do3, Prop. 4] ; we know that the dynamical system FSoHilb^g has an attractive 
fixed point hi^^^r at rank r and that the convergence of this dynamical system 
is exponentially fast. We obtain this way a new form 

= Ci{h^g^r)- 

For the second step, i.e in order to find the balanced metric /iRic-i(a;o),rj we 
introduce the operator Hilbuj^ ^, by 



e 



M 



where (po,r is the potential of the metric h^^^r S Met(L^). Iterating this proce- 
dure leads us to define at each step k a dynamical system Hilb^^ o FS which 
has an attractive fixed point huj,. ^ € Met{U') e~'^'" '^|.|o locally). Note that 
for a generic k, we have 



mib^,,. a,)) (5., 5,) = / {S^,,S,)e-'^^-^^^Lol, (44) 

Corollary 5.1. Under above assumptions, and for r sufficiently large, the se- 
quence Ci(/ii{'^J converges when r tends to infinity to the solution of the Monge- 
Ampere equation (19) with exponential speed of convergence. 
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Conjecture 1. Under above assumptions, the sequence ci{hj^^^) converges when 
I tends to infinity to a Kdhler metric LUr in the class Ci{M) with exponential 
speed of convergence. If M is Kdhler- Einstein , then uj^ converges to a Kdhler- 
Einstein metric. 

Let's describe now how our discussion can be useful for numerical approx- 
imations of Kahler- Einstein metrics on Fano manifolds. One has to notice at 
this stage that we can write (44) as 




where the (Si^k-i) G H^{M,L'^) form an orthonormal basis of holomorphic 
sections with respect to the metric Hk-i.r = Hilb^^_^ ^i^t^Jk-i r) computed 
at the previous step, i.e the Ric^^'''^^' ((jJo)"-balanced metric. Note that the 

expression (45) makes sense since the term jSi^fc-ip^ can be considered 

as a section of K^^ (g) Km^ ■ 

In order to obtain the Ric~'^([Jo)"-balanced metric Hk,r one needs to iterate 
the operator Hilbu,^^oFS, which gives a sequence {Hk^r,p)p£N G Met{H^{M, L^)) 
initialized with Hk^r.o ~ Hk~i,r- Here, we denote Hk^r = Hk,r.oc- Remark now 
that if one expects the algorithm to be convergent, and thus Hk-i,r to be close 
to Hk.r for large k and r, then it is natural to assume that Hk-i,r,i is close to 
Hk-i,r,oo = Hk,r.o, i-c just onc stcp is sufficient to get the Ric^'^'(a'o)"-balanced 
metric. This justifies at least formally the definition of canonically balanced 
metrics for Fano manifolds that appeared in [D()3, Section 2.2.2]. 

Definition 5.1. Let M be a Fano manifold and r G N* sufficiently large. For a 
metric G G Met{H"{Af -K^j), we define the operator T on Met{H"{Af -K^j) 
by 

for (Si) S H'^{A-{, —K^j) a G-orthonormal basis. A canonically balanced metric 
Gcan G M et{H^ {M , —K^j)) is a fixed point of the operator T. 

The operator T seen as acting on Bergman type metrics is actually a finite 
dimensional approximation of the Ric~^ operator. Hence, when they exist a 
priori, we expect the behaviour of the sequence of canonically balanced metrics 
to be understood via the iterations of the Ric~^ operator, and thus to be related 
to the behaviour of the Kahlcr-Ricci flow. In that direction, we can show 

Theorem 5. Assume that AI is a Fano Einstein manifold with no nontriv- 
ial holomorphic vector field. Then for r sufficiently large, there exists a se- 
quence of canonically balanced metrics Gcan,r S M et{H'^ (At , — K"^)) such that 
Ci{FS{Gcan,r)) convcrgcs when r tends to infinity to the Kdhler- Einstein metric 

^KE- 

Proof. The proof is similar to the proof of the main result of [Dol]. Actually, 
a canonically balanced metric hcan,r G AIet[~K\j) satisfies that the Bergman 
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kernel is constant on the manifold, i.e for all p e M, 



Nr 



where the sections Si E H^{AI, —K^j) are orthonormal with respect to the inner 
product 



and Nr = dim7J°(M, — A'J^). Note that one can generalize the results of Z. 
Lu [Lu] on the asymptotic of the Bergman kernel in that case. The details will 
appear in a forthcoming paper where the relationship with G.I.T stability will 



In our implementations we use this notion, i.e we consider the iterations of 
the operator T. At each step wc also choose a normalisation for the metric by 
fixing essentially its volume. Our tests on toric manifolds have shown that the 
sequence of metrics defined by (45) and (46) have similar behaviours (i.e iterat- 
ing many times Hilb^^ ^ o FS or just once) and converge to the Kahler-Einstein 
metric when it exists a priori. These procedures have the advantage to skip the 
computation of the determinant of the Fubini-Study metric as required by the 
original notion of balanced metric [Del, Luo, Zha]. Therefore they are clearly 
much more efficient and we expect that similar methods could be developed 
for Einstein non Fano manifolds. Finally, remark that even in the case of CP^ 
with the anticanonical polarisation, the sequence of balanced and canonically 
balanced metrics converge at different speeds towards the Fubini-Study metric 



5.2 The case of the projective plane blown up in three 



Let us consider the toric Fano manifold blown up in 3 (non aligned) points. 
From a result of J. Song, its a-invariant is 1 and thus it possesses a Kahler- 
Einstein metric (this is also a consequence of Tian's work of classification of 
Einstein Del Pezzo surfaces). Let us mention that the Kahler-Einstein metric 
on this manifold has been very recently studied in [DHHKW] by simulating the 
Ricci flow with PDF techniques. 

We implement our algorithms (i.e in order to find balanced and canonically 
balanced metrics) using the special symmetries of this manifold. The computa- 
tions of the points on the manifolds are relatively quick since we are essentially 
reduced to a 2-dimensional real manifold and there are some symmetries given 
by the action of Z2 (reflections) and Zg (rotations). The fan of this toric variety 
is given by the six rays spanned by 



vo = (1, 0), i;i = (1, 1), V2 = (0, 1), V3 = (-1, 0),V4 - (-1, -1), V5 - (0, -1). 



and as it is well known the polytope is actually the hexagon. There are different 
ways to choose the points on this toric manifold but we decided to just generate 
the points on one of the 6 affinc charts associated to the cone formed by pairs 
of rays (wi,t;i+i) (i.e. by defining a certain cut-off function). 




be studied. 



□ 



[Do3]. 



points 
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Our program is written in C++ (compiler gcc 3.4.6) and can be launched 
essentially with 4 different algorithms. For each algorithm we print the scalar 
curvature of the new computed metric computed at each point of the manifold. 
This gives a picture with different colours and one can easily see after some 
iterations that the scalar curvature varies essentially for the points close to the 
edges. The computation of the scalar curvature is possible with exact precision 
(up to the machine precision) since our metrics are algebraic. It has the disavan- 
tage to take time since it involves derivatives of order 4, but on the other hand 
it gives a visual output of our work, see Figure (1). Some animated pictures 
generated by the program and the program itself can be downloaded from the 
website of the author^. Despite of this loss of time, all the four algorithms (for 
the given parameters below) can be run in 1 minute or less on a decent desktop 
computer. This proves the efficiency of the methods and let us hope that it is 
possible to compute Kahler-Einstein metric on 3-folds with few symmetries. We 
now describe the results for each algorithm. 

For the computation of the balanced metric as defined in [Dol, Do2, Do3], 
we choose the parameter r = 8 (see Section 5.1) and compute approximatively 
lO"*^ points. After 50 iterations, the average scalar curvature on the manifold is 
0.95 and the maximum error is 16%. 

For the computation of the canonically balanced metric as defined in [Do3] 
or our discretization of the Ricci flow (45), we fix again r = 8 and compute 
approximatively 5 • 10** points. After 35 iterations, the average scalar curvature 
on the manifold is now 0.99 and the maximum error is less than 4%. 

Finally we try to improve the our first two algorithms by using the metric 
at rank r to compute the metric at rank r + 1. This is based on a very simple 
argument that we describe now. If one knows the balanced hr € Met{L^) (or 
canonical balanced metric) at rank r, then using the asymptotic of the Bergman 
function for higher tensor powers (see [Ln] for details) one can write on this 
surface 



where T{p) is a certain function which in fact is an algebraic expression of the 
curvature of ci{hr) and its derivatives. Here (<S'i)i=i,..,Afr is an orthonormal basis 
of H'^{M, U") with respect to the L^-metric corresponding to the choice of our 
algorithm. Once we have computed hr, it is clear that we can deduce the value 
of r at each point of the manifold. Now at rank r + 1, we look for a metric /i^+i 
such that 



with respect to the corresponding L^-mctric. Roughly speaking, it corresponds 
to force the algorithm to get a metric with constant scalar curvature up to an 

^http: / /www. ma. ic.ac.uk/ ~jkeller/ Julien-KELLER-progs.html 
Some other programs for other Fano Einstein surfaces will be available at this address. 



Nr 



scal{ci{hr)){p) 
2 



+ r(p) + 0(l/r) 
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Figure 1: Iterations no. 10,16,20,40 for the canonically balanced metric 



error of size O {l/{r + 1)^) (instead of only O (l/(r + 1)) for the balanced met- 
ric). We apply also the same trick for the canonically balanced metrics. Finally 
we call these new sequence of metrics as 1-step recursively balanced (1-s.bal in 
short) or 1-step recursively canonically balanced (1-s.c.bal in short). 

The advantage of this method is that it is particularly simple (at least for 
dimension 2) to program it, since we have already coded the computation of the 
scalar curvature. On the other hand it asks of course a little bit more time of 
computation at each step. We give now an overview of the results for r = 4, a 
choice of 10"* points on the manifolds and 15 iterations for each algorithm. 



Method (r = 4) 


Balanced 


1-s.bal 


Can. balanced 


1-s.c.bal 


Avg scalar curvature 


0.786 


0.888 


0.949 


0.984 


Max scalar curvature 


1.079 


1.041 


1.012 


1.041 


Min scalar curvature 


0.618 


0.717 


0.820 


0.867 


Time (sec) 


9.1 


16.1 


8.0 


15.5 



As an example, we obtain Figure 1 p. 23 for the canonically balanced metric with 
r = 12 and various iterations. One can see that for 40 iterations we get a metric 
with scalar curvature almost equal to 1 everywhere. Using this metric, one can 
find numerical approximations'^ of geodesies for the Kahler-Einstein metric. It 
seems to give a numerical evidence that the geodesic equations on this manifold 
form an integrable system. 

^http: / /www. ma. ic.ac.uk/ ~jkcller/ Julien-KELLER-progs.html 
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6 Iteration for other classes 



6.1 The case of negative first Chern class 

Assume now that Km > i.e that M has negative first Chern class. It is well 
known by [Ya] that such manifolds have a unique Kahler-Einstein metric ujke 
such that 

Ric{uj ke) = -^^KE- 

Let Ricw + u! = \/—\ddh with j^^e'^u)" — 1. In that case, the naive dis- 
cretization of the normalised Kahler-Ricci flow leads us to define the iteration 
scheme by 

- RicuJi = 2uji - bJi^i. (47) 
This can be rewritten in terms of the following Monge- Ampere equations 

w", ^ =e2'^'-'^'-icj",-i , ^/ e^'^'-'^'-^c^r-i = 1, VZeN. (48) 



M 



and we will denote as previously $; = Sj=i ^^^"^ uji ^ oj + \/— 109$;. Note 
that by the general results on Monge- Ampere equations [Ko], for any /, there 
exists a smooth solution Lpi of (48) and thus the iteration scheme is well-defined. 
By concavity of the log applied to (48), we have directly 

Lemma 6.1. Let /.i = —1. Then 

Eo{lji-i,i^i) = -J{LJi,uji+i) + i Jj^j{2ipi - (^z-Ow;".! < (49) 
i.e the iterations (4-8) decrease the functional Eq. 

Proposition 6.1. Let p G (l,oo). There exists a constant C\ depending only 
on p, M and lo such that 

\\H\l^[m,^)<Cu V/eN. 

Proof. Let p G (1, oo). According to the work of [Ko] (or [Bl] if p € (2, oo) which 
will suffice for the proof) and in view of (48) it suffices to prove that 

l|e''+'*'-*'-MlLnA/.c.)<C2, yieN, 

for C2 = C{M,uj,p). From the normalisation in (48) it follows that Jm(^ + 
2$; — $;_]^)u;" < where h denotes the Ricci deviation of uj. In particular 
supm(/i + 2$, - $,_i) < --i Jj^^ih + 2<Pi - - snp^ih + 2$; - <^>i-i))uj". 
Therefore 



M 



< g-PT- /j,j('i+2*,-*,-i-supj,,(h+2<I.,-$,_i))c^" 



Let G be the Green function for A = Ag satisfying J^^ G{x, y)uj^''{y) = 0. Since 
Wo + 2$; — ^i-i = 2u!i — uji-i > the Green formula gives 

^ / (sup(2$z - $i_i) - (2$i - $i^i))a;" 
V Jm m 

= -T7 f G(a;o,y)A(2$, - $j-i)c^"(y) < 
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where A = — infMxM G{x, y) and xo G M satisfies (2$;— <I>;_i)(xo) ~ supj^j(2$/ — 
Tlierefore C2 = eP(°"'=^^''+"^). □ 

Combining the previous proposition and Lemma 6.1, we obtain 

Corollary 6.1. Let M he a manifold with K]\i > 0. Then, the sequence of 
Kdhler metrics defined by (47) converges to the Kdhler- Einstein metric ujke on 
M. 

Actually, another iteration scheme is also very natural, especially if we as- 
sume the metrics uji and wi+i from (47) to be very 'close' for a large I. We can 
define the iteration scheme 

— Ric2; = (50) 
which means that we solve at each step the following Monge- Ampere equations 



h+T.-=i^,^^n^ e''+^^=i = 1, VZeN. (51) 
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Let's call = Y^j=i'fj- Up to our knowledge, it is not clear whether the 
iteration scheme (50) decreases any functional. Now, one can check that the 
proof of Proposition 6.1 can be adapted to that case without must change, i.e 
one has straightforward a C°-estimate of the potentials and thus convergence 
up to extraction of a subsequence. On another hand, the iteration scheme (50) 
is more natural from the finite dimensional perspective as we shall explain now. 

Similarly to the notion of canonically balanced metric introduced by (46) 
in Section 5.1, one can define a notion of canonically balanced for manifolds 
with negative first Chern class. Indeed, it turns out to consider the iteration 
procedure defined by the iterations of the map T = Hilb o FS where 

Hifb{h){S^,Sj)^ I /i^ «) 5, (g) S7 

for h G Met{K'^j) and Si e H'^{M, K'^j). Here we see h~k' as an element 
of Kj^ (g) Km ■ From our previous discussion p. 20 Section 5.1, one can 
expect naturally the iterates of T to approximate in the space H^{M, K\j) the 
discretization of the normalised Kahler-Ricci flow given by (50). 
As explained in [ ], another remarkable fact is that one can define an iterative 
procedure in higher dimensional spaces, i.e by defining now the maps T ~ 
rniboFS : Met{Klj) -> Met{K'll^) where 



HiTbih) {S„Sj)^ [ h® S, ® S. 

JM 



] ■ 

M 

Let's define hk = T^''\ha) for a metric /iq S Met{K^) for fco sufficiently large. 

As observed by Weinkove and Song [SW2, Theorem 1], one has a uniform con- 

1 

vergence for k +00 of the metrics to the smooth hermitian metric 

Hke G Met{Ki\j) such that ci{hKE) = ^ke is the Kahler-Einstein metric 
on M . Nevertheless this method does not seem to have numerical applications 
straightforward since it involves to consider higher rank matrices (and thus their 
inverses). 
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Remark 6.1. One can also derive an iterative scheme for manifolds with trivial 
canonical bundle. This leads to a refinement of the notion of v- balanced metric 
[Do3], and is useful for our numerical study of Ricci-flat metrics on Calabi-Yau 
threefolds. This will be addressed in a forthcoming paper. 

6.2 The case of non-canonical classes 

Let oj G Kan be a Kahler representative of an arbitrary class in the Kahler 
cone of a Fano manifold M. Let a € Ka^c^ be a representative of ci(M). By the 
Calabi-Yau theorem there exists a unique Kahler representative ui of whose 
Ricci form equals a. Define a map 

Ric[j-'^ : Kan 

by RIcq := wi (when appropriate norms are chosen, it defines again a home- 
omorphism of Banach spaces) . We may therefore define an operator Ric [^^j on 
Kan by 

Ric[w] := Ric^ ^' o Ric o Ric. 

This gives rise to an iteration on Vl by setting Ric^^ := Ric^j o Ric'-'^-' o Ric 
for each k G 1. Note that this induces again a filtration on Kan defined by 

Ka^^^ := {uj e Kan ■ Ricw G Ka^c^}. 

From what we have proved previously, we expect that the orbits of the new 
dynamical systems converge, if and only M is Kahler-Einstein, to a Kahler rep- 
resentative which is characterised by the property that its Ricci form is Kahler- 
Einstein. For Vl 7^ ci(M), one can ask how is this dynamical system related to 
the study of the space Kan- In particular, it would be interesting to relate this 
to the existence problem of extremal metrics for Kahler classes near ci{M). 
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